SYSTEM AND METHOD FOR GEOREFERENCING DIGITAL 

RASTER MAPS 

CROSS REFERENCE TO RELATED APPLICATIONS: 

This application is a continuation of and claims priority from pending 
U.S. Application "System and Method for Georeferencing Digital Raster 
Maps" (09/537,849), filed March 29, 2000. Furthermore, this application is 
related to and claims priority from the following pending applications: 
"System and Method for Performing Flood Zone Certifications" (09/537,161), 
filed March 29, 2000 and "System and Method for Synchronizing Raster And 
Vector Map Images" (09/537,162), filed March 29, 2000 which are hereby 
incorporated by reference. 

BACKGROUND OF THE INVENTION 
Technical Field: 

The present invention generally relates to geographic map processing 
and in particular to a system and method for associating digital raster maps 
with known geographic features. Still more particularly, the present invention 
relates to a system and method for associating specific points on digital raster 
maps with a geographic coordinate system. 

Description of the Related Art: 

A digital raster map is a computerized map image that resembles a 
conventional paper map in that it presents an image of the mapped area, but has 
no additional underlying data associated with the features of the map. A raster 
map is typically created by scanning a conventional paper map, and is a grid- 
based map composed of pixels (or dots) of color or black & white. Each pixel 
in the grid has can be referenced by the pixel coordinates, and has only one 
value, which indicates the color of that pixel. Raster images are commonly 
referred to as "bit mapped." 
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A vector map uses lines and polygons, rather than pixels, to describe an 
image. Unlike a raster map, in which the map image is simply stored as a grid 
of pixels, when a vector map is displayed, it is drawn from a set of underlying 
data. The vector map is created through the placements of nodes on a plane 
and connecting those nodes with lines. Vector lines can be attributed with 
tables of data such as elevations, values, names or other information relative to 
the line. Vector data can be displayed in three dimensions if the lines are 
attributed with z values, modified or changed relative to user need, or layered 
to allow for turning off and on the viewing of different information. 

Because of their feature attribution properties, they are particularly 
useful for displaying geographic data. Vector maps are used to display 
boundaries or lines that denote the position and extent of features, such as 
county boundaries or lines denoting stream and river systems. It is also very 
easy to view or manipulate the data underlying a vector map, for example to 
view or change the elevation of a feature. 

Also because vector maps are commonly drawn from a geographic data 
set/describing the area shown, they are very easily, and generally inherently, 
/eoreferenced. Georeferenclng is the process of relating source coordinates to 
referenced geographic coordinates., which are typically in standard 
latitude/longitude. An image ir a vector file is georeferenced to be used within 
a mapping/geographic environment. In a vector map, the data from which the 
map is drawn will typically already include a geographic coordinate set. 

Modern GIS systems normally make use of digital vector based map 
information. However, a vast legacy of paper-based map information exists. It 
is very expensive and time consuming to convert all of the information on these 
paper maps over to a digital vector format. In many cases the scope and 
expense of such conversions render them completely impractical. However, 
even when a complete conversion to digital vector-based format is not possible, 
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it is still possible to obtain some of the benefits of computerized map systems, 
first by converting the paper maps to digital raster maps (by scanning them), 
and then by georeferencing the raster image. After georeferencing, there should 
be a clear relationship between the pixel coordinates in the raster map, and the 
5 geographic coordinates of the feature represented by that pixel. It would be 

desirable to provide a system and method for georeferencing a raster map by 
associating points on that map with corresponding points on a previously- 
georeferenced vector map, or with reference points designated in the raster map 
which have known latitude and longitude. 
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SUMMARY OF THE INVENTION 

It is therefore one object of the present invention to provide improved 
geographic map processing. 

It is another object of the present invention to provide an improved 
5 system and method for associating digital raster maps with known geographic 

features. 

It is yet another object of the present invention to provide an improved 
system and method for associating specific points on digital raster maps with a 
geographic coordinate system. 

10 The foregoing objects are achieved as is now described. The preferred 

embodiment provides a system and method for georeferencing digital raster 
maps by associating specific points on the raster map with corresponding 
points on a georeferenced vector map or another georeferenced raster map. 
According to the preferred embodiment, a raster map and a corresponding 

15 vector map are simultaneously displayed to a user. The user then locates a 

common geographic point or feature on each map, and marks each of them as a 
unique point-pair. When the user has marked at least two point-pairs, the 
system then computes a georeferencing function, based on the pixel- 
coordinates of the points marked on the raster map and the corresponding 

20 geographic coordinates of the points on the vector map. Thereafter the 

geographic coordinates of any point on the map may be easily computed. The 
preferred embodiment provides that as more point-pairs are identified, the 
georeferencing function is modified for increased accuracy. 

The above as well as additional objectives, features, and advantages of 
25 the present invention will become apparent in the following detailed written 

description. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features believed characteristic of the invention are set forth 
5 in the appended claims. The invention itself however, as well as a preferred 

mode of use, further objects and advantages thereof, will best be understood by 
reference to the following detailed description of an illustrative embodiment 
when read in conjunction with the accompanying drawings, wherein: 

Figure 1 depicts a data processing system in accordance with a 
10 preferred embodiment of the present invention; 

Figure 2 is an image of an exemplary raster map, in accordance with 
the preferred embodiment; 

Figure 3 is an image of an exemplary vector map, corresponding to the 
raster map of Figure 2, in accordance with a preferred embodiment of the 
15 present invention; and 

Figure 4 is a flowchart of a process in accordance with a preferred 
embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 



With reference now to the figures, and in particular with reference to 
Figure 1, a block diagram of a data processing system in which a preferred 
embodiment of the present invention may be implemented is depicted. Data 
processing system 100 includes processor 102 and associated L2 Cache 104, 
which in the exemplary embodiment is connected in turn to a system bus 106. 
System memory 108 is connected to system bus 106, and may be read from and 
written to by processor 102. 

Also connected to system bus 106 is 1/0 bus bridge 110. In the 
exemplary embodiment, data processing system 100 includes graphics adapter 
118 connected to bus 106, receiving user interface information for display 120. 
Peripheral devices such as nonvolatile storage 114, which may be a hard disk 
drive, and keyboard/pointing device 116, which may include a conventional 
mouse, a trackball, or the like, are connected to 1/0 bus 112. 

The exemplary embodiment shown in Figure 1 is provided solely for 
the purposes of explaining the invention and those skilled in the art will 
recognize that numerous variations are possible, both in form and function. For 
instance, data processing system 100 might also include a compact disk read 
only memory (CD-ROM) or digital video disk (DVD) drive, a sound card and 
audio speakers, and numerous other optional components. All such variations 
are believed to be within the spirit and scope of the present invention. Data 
processing system 100 is provided solely as an example for the purposes of 
explanation and is not intended to imply architectural limitations. 

The preferred embodiment provides a system and method for 
georeferencing digital raster maps by associating specific points on the raster 
map with corresponding points on a georeferenced vector map or another 
georeferenced raster map. According to the preferred embodiment, a raster 
map and a corresponding vector map are simultaneously displayed to a user. 
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The user then locates a common geographic point or feature on each map, and 
marks each of them as a unique point pair. When the user has marked at least 
two point-pairs, the system then computes a georeferencing function, based on 
the pixel-coordinates of the points marked on the raster map and the 
5 corresponding geographic coordinates of the points on the vector map. 

Thereafter the geographic coordinates of any point on the map may be easily 
computed. The preferred embodiment provides that as more point-pairs are 
identified, the georeferencing function is modified for increased accuracy. 

Figure 2 is an exemplary raster map, in accordance with the preferred 
embodiment. This exemplary map shows a scanned image from a Federal 
Emergency Management Agency (FEMA) paper map. This raster image shows 
a land area with flood zone indications, but would, in a computer system, 
contain no underlying data regarding the area shown. 

Figure 3 is an exemplary vector map, corresponding to the raster map 
of Figure 2, in accordance with a preferred embodiment of the present 
invention. This map shows the same area as the map in Figure 2, but is created 
by a computer system from a database describing the locations of features such 
as the streets shown. Typically, each feature shown on a vector map such as 
this will already be georeferenced, in that the geographic coordinates of each 
feature will also be recorded in the underlying data 

A digital map image is considered georeferenced if a pair of 
mathematical functions, f, and g, have been defined that can be used to convert 
back and forth between the coordinates of the map image (as defined by the 
pixels of the image) and the corresponding longitude and latitude of the 
25 location of that point. That is, f and g do the following: 

1 . If (x, y) represents a location on the digital map image, then f(x, 
y) =(Lon, Lat) represents the longitude and latitude of the corresponding 
physical location. 
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2. If (Lon, Lat) represents a physical location that lies within the 
region covered by the map, then g (Lon, Lat) = (x, y) represents the point on 
the digital map image that corresponds to that longitude and latitude. 

Here, x and y represent the natural internal coordinate system of the 
5 map image. Typically, as described above, a digital raster map image uses the 

pixels of its image as a natural coordinate matrix. However, in most cases, a 
vector-based map image uses longitude and latitude as its internal coordinate 
system; if so, it can be considered to be trivially georeferenced already. 
Therefore, the functions f( ) and g ( ) above are non-trivial 

10 georeferencing functions required to convert back and forth between coordinate 

systems. 

Following is a description of a georeferencing process in accordance 
with the preferred embodiment: 

First, the raster map to be georeferenced (referred to as Mapl) is 
15 displayed on a computer monitor in conjunction with another previously 

georeferenced map (referred to as Map2) (step 400). Typically, Map2 will be a 
vector map, but could also be another raster map, for which georeferencing 
functions are already known. Mapl and Map2 can be independently 
manipulated for rotation, position and scale. The user will manipulate the maps 
20 until the region shown in Mapl is contained within the region shown by Map2. 

Next, pairs of points representing identical geographic locations on 
Mapl and Map2 are identified, and the user uses a keyboard, mouse, or other 
pointing device to mark these point-pairs (step 405). For each point-pair, one 
point is marked on Mapl and the other point is marked at the corresponding 
25 location on Map2. Examples of common georeferencing point-pairs include 

street intersections, places where rivers cross roads, mountain peaks, building 
locations, and other readily identifiable map features. If any points on Mapl 
have points of known longitude and latitude already marked on the map, then a 
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georeferencing point-pair can be obtained by marking that point on Mapl and 
moving to the point with identical latitude and longitude on Map2 and marking 
it as the other half of the point-pair. Alternatively, once the known points are 
marked with the known latitude and longitude, they can be used to determine 
5 the georeferencing functions without requiring a corresponding mark on Map2. 

After a point-pair is marked, the x-y point on Mapl is assigned the 
corresponding latitude and longitude of its matching point on Map2 (step 410). 

When two georeferencing point-pairs are marked, the system will 
compute georeference functions for Mapl based on a linear transformation that 
10 allows an arbitrary rotation and predefined scaling, between the x-y domain 

and the latitude- longitude domain (step 415). 

To improve accuracy and allow increased generality, more 
georeferencing point-pairs will frequently be desired. To facilitate the process 
of gathering these additional georeferencing point-pairs, Mapl and Map2 are 

15 now selectively synchronized. This means that as either map is zoomed, 

panned, scrolled, or otherwise caused to display a different region, then the 
other map automatically does the same. This greatly increases the ease with 
which georeferencing point-pairs can be identified and marked on the maps, 
since when the user of the system locates an identifiable feature on one map, 

20 the other map will automatically be showing very nearly the same location on 

the other map. 

Thereafter, when the user marks a georeferencing point on one of the 
maps, then the system automatically marks the corresponding point on the 
other map, based on its current georeferencing function (step 420). The user 
25 can then make corrections to the location of a point marked by the system in 

this way (step 425). The user may enter as many point-pairs as are desired; 
each additional point-pair will increase the accuracy of the georeferencing 
process. 
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When three or more georeferencing point-pairs are determined, the 
system computes georeference functions for Mapl based on a completely 
general linear transformation between the x-y domain and the latitude- 
longitude domain. 

When four oV more georeferencing point-pairs are determined, the 
neral linear georefeiencing functions are over-determined. This means that 
more than the required amount of information to compute the general linear 
georeferencing functions is available, but that it is not, in general, completely 
consistent. The system ike the extra information contained in the additional 
georeferencing points to\ provide validation checks to protect against the 
possibility that some of tHp data points may be inaccurate (step 430). Points 
that deviate excessively with respect to a calculated standard error are 
presumed to be inaccurata and are omitted from the calculation of the 
georeferencing functions. Noie that as new point are added, the system also 
rechecks points previously maMced as inconsistent, to determine if those points 
should now be considered wheiMecomputing the georeferencing functions. 

The additional point-pairs allow the system to compute the general 
linear georeferencing functions which best fit the combination of all the 
available georeferencing point-pairs, so that the active georeferencing functions 
are revised with each new point-pair (step 435). According to the preferred 
embodiment, this is done by using a "least square" parameter fitting operation. 



The userWay then proceed o enter the next point-pair (step 440) . 
yWhen the user is finished, the system stores the active georeferencing functions 
/ with the raster-mapVstep is 445). At this time, the raster map is considered 
/ fully georeferenced. When accessed at any future time, the system may simply 
retrieve the georeferencmg functions, and apply them to find the latitude and 
longitude of any point on the raster map. 

The process of determining he georeferencing function set from a set of 



108344.00016: 2613904 



10 



10 




20 



point-pairs is believed tAbe within the ability of one of ordinary skill in the art. 
The specific approach msed by the system an method of the preferred 
embodiment is discussed below. 

At any given point in the georeferencing process a set of points P, has 
been specified. Let / e P denote one of the points of P . Associated value x i , 

y g9 Lat t and Lon i refer to thejc and y coordinates of the point on the bitmap, 
and the latitude and longitude of the point on the earth's surface. We also 
define AcPto be the subset of "active" points, i.e. points which are 
currently being used in the calculation of the georeferencing functions. Assume 
that we have selected a parametric family of functions F. Functions in this 
family map from (x 9 y) to(Lon,Lat). Any function / e calFhas components 
f{x 9 y) = Lon 9 and f 2 (x 9 y) = Lat . From within this family we seek the 

specific function / , which (in a sense to be defined below) comes closest to 
making the following system of equations true: 



/(*/ >yt ) = {Lon^h^ ) for i e A 



nee determined, / will b the georeferencing function which is used 
pute corresponding latitude and longitude values, (Lon 9 Lat) for any 
{x 9 y) on th^ bitmap. There are any number of possible ways to define 
the function that "cbmes closest to making (1) true." We shall follow a "least 
squares" approach \so known in mathematics as an L 2 approach. This 

approach seeks to fin)^ the function, /, which minimizes the sum of the 
squared differences betW^en the actual and the predicted values of latitude and 
longitude. In other words,^rom among all the functions f e F , /is the one 
which minimizes: 



25 



(2) 
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brrors, or so that it m: 



10 




Among varibus alternative methods for choosing the function / are 
choosing it so t ht it minimizes the sum of absolute errors (rather than squared 

to* 



es the largest error. Other criteria are also possible. 




In what follows we shall describe methods for computing 
georeferencing functions based on both a "rotational linear" parametric family 
of functions and a "general linear" family of functions. A knowledgeable 
practitioner will readily perceive that this approach may be easily generalized 
to other parametric families of functions. Obvious candidates include nonlinear 
parametric families of functions arising from the projection of the earth's 
surface onto a flat map. 

The General Linear Case: In the general linear case, we let F be the 
of all possible linear transformations which map from {x 9 y) to (Lon,Lat). 
Thus 







X 












+ 




a 22. 









(3) 



15 for some choicAof parameters d u 9 d 129 d 219 d 22 9 £ 19 and b 2 . If the 

region covered by the mab to be georeferenced is not too large, then this family 
of functions will contain k suitable function , / whose total error is quite 
small. In the case where the map to be georeferenced covers a larger area than 
this, then the curvature of the \arth must be taken into account and F is not a 

20 suitable family of functions. In shch a case, nonlinear functions must be used as 

mentioned above. We shall no\ pursue that case further, since it is a 
straightforward extension of the procedures used in the linear case. 




To find / we seek me parameters which minimize 
SSE = X(fl„x, +a n J'\*/ -Lon.f +(a 21 x f + a 22 y, +6 2 -Lat,) 2 The 
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parameter values which minimizk this expression are found by solving the 
following two independent systems of linear equations: 
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These systems can be easily solve by well known methods, such as 
Jaussian Elimination, or LU factorization. The solutions yield the desired 

values ofa n , a n , a\ a 22 , b x , and b 2 . It should be noted that equations (5a) 
and (5b) do not have a Wque solution unless three or more non-co linear points 
are contained in A . Generally speaking, then, it requires 3 points to choose a 
georeferencing function from the family of general linear transformations. 
When there are four points or more, it is possible to compute a standard 
deviation of errors using the formula: 



s = 



Z[( 5 n*/ +5 i2X +X-£o7i,.) 2 +{a 2l Xi +b 2 -Latj] 

reA \ 



V«-3 



(6) 
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s is an estimator for the amount of error to be expected between actual 
and predicted latitude and longitudes values. 

Note that the inverse georeferencing function, / that maps from 

(Lon, Lat) to (x,y) is readily obtained now, by inverting the function /. 
Having done this, it is possible to compute a similar standard deviation of error 
for f~ l which is an estimate of error as measured in the bitmap coordinates. 

The Rotational Linear Case: In the, so called, rotational linear case, 
we let Fbe the set of all possible linear transformations which map from 
(x,y)to {Lon, Lat) and which also allow for a known scale change caused by 
the difference in distance per degree of latitude and per degree of longitude at 
any point not on the equator. Thus . 



'A. -a' 


X 




"A" 
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-A -A. 









(7) 



15 



for some 

choice of parameters fi x , y£ 2 , y# 3 , and >ff 4 , where y =cos(Lat) is 

computed from the estimated latitude of the map to be georeferenced. 
Clearly, this is a special case of the more general linear transform 
described above, where we take 

a n =& / r>a X2 = -fl 4 > a 2i ="A> 



20 



In theory, if the region covered by the map to be georeferenced is not 
too large, then this family of functions will contain a suitable function, 
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/ whose total error is quite small. In the case where the map to be 
georeferenced covers a larger area than this, then the curvature of the 
earth must be taken into account and F is not a suitable family of 
functions. In such a case, nonlinear functions must be used as noted 
above. In practice, the general linear transform will yield results which 
are somewhat superior to the rotational linear transform due to the 
former's innate ability to use arbitrary scaling factors. These arbitrary 
scaling factors might compensate, for example, for scaling flaws in the 
digital scanner used to scan in the maps. The primary advantage of the 
rotational linear transform is that it can be computed with only two 
reference points (unlike the general linear transform which requires 
three reference points). 

To find / we seek, as before, the parameters which minimize SSE as 
defined in equation (4). 
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These systems can\easily solved by well known methods, such as 
Gaussian Elimination, or LO factorization. The solutions yield the desire 

values of fi x , J3 19 ft X9 and y9\which in turn yield the desired values for a,,, 

<*\2> <*2\> <*22> *1 > ^ b 2 . 



10 



It should be noted that equation (8) does not have a unique solution 
unless two or m re points are contained in A. Generally speaking, then it 
requires two points to determine a georeferencing function from the family of 
rotational linear transformations. WhenUhere are three points or more, it is 
possible to compute a standard deviation of error, s using the formula: 



^[(flnX, +5 12 J>, +b x -LonJ +(fl 21 x, +a 22 y t +£ 2 -Lat J 



s = 



1 leA 



n-2 



s is an estimator for the amount of error to be expected between actual 
and predicted latitude and longitude values. 

Note that the inverse georeferencing function f~ l that maps from 

15 (Lon, Lat) to (x,y) is readily obtained now, by inverting the function /. 

Having done this, it is possible to compute a similar standard deviation of error 
for f~ l which is an estimate of error as measured in the bitmap coordinates. 

Automatic Error Detection and Handling When individual points 
are being assigned k y, Lon, and Lat values, there is always a potential for 
20 ^)>S y^rror. To reduce the rkk of incorrect georeferencing resulting from such errors, 
certain error handling procedures are built into the georeferencing process. The 
fundamental conk t is tn^t of detecting a "bad" point and then removing it from 
the set of active points, A Note that removing point from A will not delete the 
information associate, witji that point, but it will cause the georeferencing 
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parameters to be completely uninfluenced by that point. We o not wish to 
remove the point entirely, smce it may be determined at a later stage of the 
georeferencing, that the point was not really bad at all, and should be-used in 
the georeferencing calculation. This will be clarified shortly. 



"^S. Detecting Ikid Points The following steps outline the bad point 
defection process \using the general linear transform approach to 
/georeferencing. 

1. Begin by placing all existing points into the active set, A 

2. If there are fewer than five active points then you are done, 
Otherwise, for each oVthe currently active points in turn, move it (call 

it point k for the sake of convenience) temporarily out of the active set, and 
then calculate the resulting inverse georeferencing function (call it g^) and 
its corresponding SSE k . AIsol calculate the difference between the 

predicted value and the actual\ value S k = \g^(Lori k , Lat k ) - (x k , y k )| . 

Make a note of the values, S k and • Return point k to the active set 

(and move on to the next value of k . 

3. From among the results found in step 2 above, find the 
point, k, with the largest value of ^y^gg )c x andS k which also 

S / 

satisfies )c l andS k )c 2 , where c x , andc 2 are some constants which are 

set according to the general level of accuracy to be expected on the particular 
type of map which is being georeferenced, the current number of active points, 
and the dots per inch of the scanned image. If there is such a point then mark it 
as bad (by removing it from the active set) and return to step 2 above. 
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Otherwise you are done. 



There ara several things to note about this procedure. One is that it 
Jlowing the value of c x and c 2 to change with the number of active points, 
makes it possible i©r the georeferencing system and method to utilize points 
which it might originally determine bad or inconsistent after a large enough 
sample of points hasXbeen gathered to make it clear that a lesser level of 
accuracy is all that cam be achieved on this map. Another observation is that 
using this procedure it re impossible to reduce the number of active points 
down to less than four (unless you started with less than 4 in which case this 
procedure does not apply at Vll). This scheme means that as each new point is 
added, all points determines so far are considered, even those which had 
previously been marked bad. Vhus early "misjudgments" on the part of the 
system can be corrected later, in PUjht of new point information. 

The same bad point detection process, can also be implemented using 
the rotational linear transform approach. In this case the method is capable of 
reducing the number of active points down as low as three (rather than four for 
the general linear transform approach outlined above). This can be useful when 
dealing with small sets of active points. 

When the system is georeferencing raster map images that cover a large 
enough area so that the (nonlinear) curvature of the earth is a source of 
significant error, the system can calculate nonlinear georeferencing functions 
whose form corresponds to the map projection that was used to create Mapl. 

A specifics example of the operation and application of the preferred 
georeferencing method may be shown with reference to the "Flood Zone 
Determination" business. The Federal Emergency Management Agency 
(FEMA) publishes a Mibrary of tens of thousands of paper maps showing 
various types of flood axmes and their locations in the United States. A flood 
zone determination on a property is frequently done in the following way: 
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1 . The address of the property is examined, and the 
location of thd property is determined (perhaps through the 
use of a geocoaing system, or b examining an available 
street map). \ 

5 2. A map analyst attempts to determine which of the 

many thousands o UFEMA flood maps will contain this 
property. \ 

3. The map analjtet goes to a map storage area and 
retrieves the desired majp, often examining several maps 

10 before making a final selection. 

4. Having retrieved thVpaper map, the map analyst 
next determines where, precisely, the property is located on 
the map. \ 

5. Finally, the map analyst Examines flood zone 

15 notations on the map at the propeiw's location in order to 

determine its flood-zone status. \ 

When performed using paper maps, the above process is difficult and 
quite time consuming . A database of scanned raster map images, alone, can be 
20 used to reduce the time and effort expended on step 3 above. However, 
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georeferenced raster map images can also be used to eliminate virtually all 
work from steps 2 and 4 above. Since those are generally the most time 
consuming steps under the current process, the value of georeferenced flood 
maps is considerable. 



Using georeferenced flood map raster images, steps 2 and 4 above, are 
^ ft /"eplaced by: \ 

%^ ( 2. A computer system combines the pre-designated outlines of 

the raster mip and the georeferencing information to obtain a polygon 
expressed in Wins of latitude and longitude that outlines the 
10 region includes in each flood map. Then the system determines which 



of the polygons\contain the address in question, which is done using 
a "point-in-polygpn" algorithm. At the conclusion of this process, the 
computer system nas identified a map panel (or perhaps a small 
number of map panels) that contains the address. 

15 4. Since the latitude and longitude of the property are known (by 

virtue of a geocoding\phase), the computer system can use the 
georeferencing of the map panels to locate the property on each of the 
panels found above, tluA largely eliminating any need for he map 
analyst to scan the flood map for the address location 

20 Therefore, it is easily seen that by using georeferenced raster maps 

instead of paper maps, the process of determining flood zones is substantially 
automated and much more efficient. 

It is important to note that while the present invention has been 
described in the context of a fully functional data processing system and/or 
25 network, those skilled in the art will appreciate that the mechanism of the 

present invention is capable of being distributed in the form of a computer 
usable medium of instructions in a variety of forms, and that the present 
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invention applies equally regardless of the particular type of signal bearing 
medium used to actually carry out the distribution. Examples of computer 
usable mediums include: nonvolatile, hard-coded type mediums such as read 
only memories (ROMs) or erasable, electrically programmable read only 
5 memories (EEPROMs), recordable type mediums such as floppy disks, hard 

disk drives and CD-ROMs, an transmission type mediums such as digital and 
analog communication links. 

While the invention has been particularly shown and described with 
reference to a preferred embodiment, it will be understood by those skilled in 
10 the art that various changes in form and detail may be made therein without 

departing from the spirit and scope of the invention. 
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